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Abstract 

We study the dynamical behaviour of the proton transfer in the hydrogen 
bonds in the base-pairs of the double helices of the DNA type. Under the 
assumption that the elastic and the tunnelling degrees of freedom may be 
coupled, we derive a non-linear and non-local Schrodinger equation (SNLNL) 
that describes the concerted motion of the proton tunnelling. Rough esti- 
mates of the solutions to the SNLNL show an intimate interplay between the 

concerted tunnelling of protons and the symmetry of double helix. 
PACS numbers: 87.15-v, 87.14.Gg 
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I. 

The recent direct observation of coherent proton tunnelling in macromolecules ]I] has 
focused the attention on systems allowing for the nontrivial dynamics of protons contained in 
hydrogen bonds. Perhaps, the most significant system of the kind is the DNA molecule. The 
problem of the proton transfer inside the hydrogen bonds of the double helix representing 
a molecule of DNA, has a long and rich history. The importance of the phenomenon was 
noticed soon after the DNA double helix had been discovered (see ||). In fact from the 
chemical point of view the proton transfer is a so-called tautomerization reaction, that is a 
kind of transition that preserves the constituent atoms of a compound but at the same time 
changes their mutual positions. It has been believed that the tautomerism transition could 
provide a mechanism for genetic mutations ||. 

Later Crick |J suggested that the mutations could be due to conformational changes 
within the double helix, the so-called "wobbling" ||. By now there has been no definite 
conclusion as to whether the wobbling or the tautomerism are responsible for the occurrence 
of mutations. The recent work on the i-motif of DNA, in which the tautomerism, rare in the 
usual DNA, is the rule, and not an exception || has given a new impetus to the problem of 
proton transfer. 

The phenomenon of tautomerism in the base pairs of DNA has been studied extensively, 
both experimentally and theoretically @, || , |J, but for individual single molecules of 
purines and pyrimidines. Here it should be noted that even though one has got a considerable 
body of information about the chemical reaction corresponding to the tautomer transition, 
such as its constant of reaction, and even managed to find an estimate for the concentration 
of tautomeric forms, 10~ 4 — 10~ 5 mol/1, for the B-type of DNA, there is a little knowledge 
as far as the dynamics of the proton transfer accompanying it in DNA is concerned. In fact, 
the estimates for the transition frequency widely diverge, and generally it is believed to be 
within 10 6 - 10 11 a" 1 , §. 

In our opinion the unsatisfactory state of the art as regards to the proton transfer, and 
the tautomeric dynamics, in the DNA, is to a certain extent due to the absence of theoretical 
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models of the phenomenon. In the present paper we would like to draw attention to the 
fact that if the tautomeric transitions are coupled with a change in elastic properties of 
the DNA molecule, one could expect a concerted tunnelling of the protons in the hydrogen 
bonds. Taking into account the extremely sophisticated nature of the system, i.e. the DNA 
molecule, we aim at studying it in a qualitative framework, which takes into account, the 
main features of the molecule, i.e. the presence of the two strands, the helical structure, or 
the winding symmetry, and the dynamics of protons considered within the approximation of 
two-level systems. 
II. 

Aiming at a simple model, we shall describe the states of the hydrogen bonds as those 
of a one-dimensional Bose oscillator, the two-level requirement being accommodated by 
considering only its ground state and the lowest excited one. In this sense the tautomeric 
reaction of the proton transfer is described by the excited state of an oscillator. Thus, we 
shall consider the protons inside the hydrogen bonds of a molecule of DNA as a quantum 
system described by the Bose operators b+,b n subject to the usual commutation relation 

in which n, m are the indices of the corresponding sites of hydrogen bonds. We shall suppose 
that the system of the protons is in a weakly excited state that can be cast in the form of 
wave function 

|£>(f)>=E4,(f)tf|0> (1) 

n 

in which A n (t) are complex amplitudes subject to the relation 

EKI 2 = i 

n 

The Hamiltonian of the total system, that is the protons and the elastic part corresponding 
to the sugar-phosphates, reads 

H = H H + H Y + H I (2) 
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where Hh is the proton Hamiltonian 



h h =y: E b + n b n - K ^(b+ +1 b n + b+b n+1 ) , 



(3) 



n n 



where Eq is the level splitting between two states of the proton in a hydrogen bond, and k is 
the tunnelling probability. Hy is that of the elastic part, and Hj is the interaction between 
them. We shall consider the elastic part of the system as the classical one, and even neglect 
its kinetic energy. The reason is that the characteristic times for the elastic vibrations of the 
DNA molecule are usually estimated as being in the region of 10~ n — 10 _13 s, whereas the 
proton tunnelling is alleged to be within 10~ 6 — 10~ n s (more exact figures are unavailable; 
according to || the value is 10~ 6 s). Thus, one may suggest that the elastic motion is 
following the proton tunnelling without inertia so that we may take into account the lattice 
deformation only through the potential energy Hy, neglecting the kinetic one. In writing Hy 
we shall follow the method worked out in papers 0,0, aiming at a simplified description of 
the dynamics of the double helix which is considered as a one-dimensional lattice of vectors 
y n that describe the mutual positions of the two strands at sites corresponding to the base 
pair of index n. The helix structure is described with the help of the covariant derivative for 
the description of deformations resulting from the change in the positions of the strands, so 
that the potential Hy reads 



Here the harmonic term (l/2)ey^ describes the binding of the two strands, and Vy n is the 
covariant derivative in discrete form, linearized with respect to rotations that relates the 
two adjacent sites, n and n + 1. 




\K 



(4) 




1 



(5) 



and 





(6) 



4 



where = (0 , , Q), K , e are elastic constants, and a is the spacing of the lattice of the 
base pairs. We are assuming that the molecule is parallel to the z axis. Having chosen the 
dynamical variables 6+ ,b n , y n for the protons and the elastic excitations, we have a limited 
number of option for the interaction energy. The simplest one should couple the proton 
excited states given by b+b n , and the elastic deformations, which are a function of y n , in 
the simplest situation we adopt, the linear one. We assume that y n enters in the form of 
its covariant derivative, for this way one might accommodate the interaction of the proton 
tunnelling with the stack of the ir— electrons of the bases. 

The choice of the interaction Hamiltonian is delicate but the knowledge of the structure 
and binding in DNA gives some indications. The stacking interaction between adjacent base 
pairs is strongly affected by the overlap of the ir electrons of the bases. But the proton 
transfer is accompanied by a redistribution of the electrons on the bases so that the stacking 
interaction is changed and therefore we expect that the proton tunneling will affect the 
coupling between adjacent bases. It means that the derivative Vy n rather than y„ itself 
should enter in the interaction Hamiltonian. We have chosen the simplest form coupling 
the commponent of the gradient which is along the direction h n of the hydrogen bonds 
connecting the bases within pair n and the state of the protons in this pair given by b^b n , 
i.e. 

#/ = -A]T(Vy„-h n )&+5 n , (7) 

n 

in which the vectors h n are subject to the helical symmetry of the molecule. They can be 
written in the form 

h n = (cos na ; sin na ; 0) , (8) 

where the angle a is related to the helical pitch. Having chosen a primitive model of the 
DNA, and neglecting the subtle features like 2 or 3 hydrogen bonds for the A -T or G - 
C base pairs, and confining ourselves to one-site one H-bond picture, we may set a = Q 
equalizing the pitch of the h n and the twist of the double helix. 
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It is important to notice that the sign of A will be crucial for the SNLNL equation that 
we shall derive. Our choice means that the stacking energy is reduced when protons are in 
the excited state. 

We employ the method worked out by Davydov [[| to study the model described above. 
According to 0, 

• we have to calculate the effective potential 

U eff =< D(t)\H Y + H T \D(t) > (9) 

which is a function of y n , 

• find its minimum, y°, 

• substitute y° into H given by Eq. (@), thus obt aining the effective Davydov Hamilto- 
nian Hp. 

Finally, we have to write down the Schrodinger equation for the function \D(t) >, given 
by Eq. ([!]), and the Hamiltonian Ho 

ih^\D(t) = H D \D(t) > (10) 

Both sides of the equation indicated above are linear forms in the operators &+, and by 
equating the coefficients at corresponding &+, we obtain the following equation for A n 

ih-^ = E A n - K(A n+1 + 4,-1) - -\A n \ 2 A n \A m \ 4 ]A n + (11) 

E cosl m '- m "l0cos[(m'-m / ')(0-fi)]|A m ,| 2 |A m ,fMn + 

ml ,m" 

^^{E cosl—l 0cos[(m - n)(0 - n)]\A m \ 2 }A n , 

in which the angle is determined by the equation tan0 = Q and X/K and ea 2 / (KQ 2 ) 1 
are small parameters. 
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In assessing the importance of cos- terms in the equation given above, it is worthwhile 
to note that since we have assumed a = Q and there is the relation tan</> = Q, the typical 



term in the equation ( LTD contains a factor that reads asymptotically 

exp(— Q 2 \m — n\)cos(\m — n\Q 3 ) 
We have assumed Q being a small parameter, and therefore the oscillations due to Q 3 are 



negligible. Let us consider a possible simplification of the equation (|TT|) obtained above. 
There is a chance that it may have a bearing on the dynamics of the proton transfer in 
DNA. To that end note that the angles <j) and Q are close to each other; indeed, for the B 
DNA the angle Q corresponds to the pitch, i.e. 10 steps for 27r. For this reason, we set all 
the functions cos[(m — n)(</> — Q)} equal to 1 in (pTTj) , and introduce 

A(z,t)=eM-i I ^)B(z,t) (12) 

Then, in the continues notations, (i.e. for scales larger than a) we obtain the following 
non- linear and non-local Schrodinger equation (z = x/a). 

.dB d 2 B A A 



dt dz 2 



^lu t \B 2 \B - -co T B j \B(z')\ A dz' + (13) 
K K J 



A ea 

-UJrp- 
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K KQ? 

A 



[I dz' j dz" exp(-fi\z' - z"\)\B(z')\ 2 \B(z"\ 2 )B(z) + 



K UT lk? { I dz '^-^ z ~ A)\B{z'?)B{z) , 



in which 



1 VL 2 
H= -lncos0 = -ln(l + fi 2 ) ~ — , (14) 

and ujh = k/H, ujt = X/h. We may perform a very crude estimate so as to see the part 
played by the non-local terms as regards the structure of solitons that might turn around. 

The standard procedure to treat the non-linear Schrodinger equation is as follows ||, 
Til . We are looking for the solution in the propagating wave (soliton) form 



B = exp[i(kz — ut)]if)(z — vt) (15) 

Introducing ( |T5| ) into ([U|) we get two equations for the imaginary and reel parts of the 
solution fll^p . From the former one we find the velocity of the soliton v 

v = 2u H k, (16) 

and the real part of fll3l) leads to the equation 

^ = -u H [-k 2 +],]- _ J \^ dz ' + (17) 

Jc<*rjjfp4> J dz 'J dz''eM-^W-z''M(z')\ 2 mz'')\ 2 + 
A f 

—u T 4> J dz'exp(-n\z - z'\)\i/j(z')\ 2 

Neglecting \x we get the standard non-linear Schrodinger equation however with renormalized 
coefficients. It is easy to find the first integral of the equation, i.e. the energy W 

W 2 + l-k 2 + v - ^(-^C + TM 2 + = w (lg) 

Where we used the following notations: 

£ = ^=-,r=/^W)| 4 

and C is a factor of the order of 1. It is worth noting that even though we have made 
rough simplifications, as regards the non-locality, its bearing on the proton dynamics still 
has remained, as is seen in the ea 2 /KQ 2 term, preserved in the equation given above. One 
may infer from the fact that under approximations used, there is a profound interplay of the 
dynamics of proton transfer and the conformational structure of the double helix. 

From the equation ( [P7D we obtain the asymptotic width A of the soliton which reads 



ujj. ea 2 



-1/2 

(19) 
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The frequency uo ac = K/h is generally accepted to be within 10 11 — 10 13 Hz, but the available 
estimates for u T ,u H ,e widely diverge, u H = 10 6 — 10 11 Hz, uj t /ojh = 0.1 — 0.01, and 
ea 2 /KQ 2 ~ 0.1. Consequently, the width of the solution varies within 10 — 1000A, so that 
one may expect a concerted tunnelling of protons for the lower estimation of A. 
III. 

Concluding, we should like to draw attention to the fact that the existence of an appre- 
ciable interaction between the proton transfer inside the hydrogen bonds of the double helix, 
and elastic modes of the latter, could result in a concerted dynamics of the protons which is 
generally of non-linear character and governed by the non-linear and non-local Schrodinger 
equation. To our knowledge, it is for the first time one is running across the non-linear 
Schrodinger equation with non-local terms. The transfer of protons may be due to various 
reasons, and among these are the action of external agents, especially, according to recent 



experimental work |TT| , |T2| enzymes; it is also worth keeping in mind the possible links 
with the mutation mechanism, In this paper we have tried to put these ideas in more 
quantitative form within the framework of a model that accommodates the basic symmetry 
structure of the double helix, and as was shown above, allows for certain rough estimates 
of dynamical features that may surface; presumably of the soliton nature. In this respect, 
it is worthwhile to note that even the crude estimate we made, conserves the bearing of 
the double helix, as in seen through the occurrence of the characteristic parameter e/ (KQ 2 ) 
in the final formulas. Of course the equation ( |TTD we obtained has a larger scope. It may 
not imply the existence of solitons in DNA, but could open another possibility, the possible 
existence of nonlinear localization leading to a collapse of an initially broad excitation into a 
higly localized deformation fll3| . Collapse does not occur in the standard NLS equation, but 



the existence of nonlinearities with higher power in our NLNLS suggest that it could occur 
in this equation, although, as these terms are non local, no definite conclusion can be given 
without further studies of the equation. If this hypothesis would be confirmed, a weak and 
broad perturbation of the hydrogen bonds of a DNA molecule by the vicinity of an enzyme 
carrying local charges could trigger this nonlinear localization phenomenon and finally lead 
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to the formation of a tutomerized form by the tunneling of one proton. As the present stage 
of the study, this is however only a speculation raised by the form of the NLNLS equation 
that we derived. 
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